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The calculation of control parameters for a system control method is based on 
the model of the system with assumed fixed internal parameters. However, 
these parameters can vary greatly due to several phenomena. This paper 
presents an adapted control of a doubly fed induction generator machine 
robust against the rotor resistance variations of the machine used as a 
generator in wind energy conversion systems. The adaptation is ensured by a 
system allowing to identify in real time the value of the resistance, the system 
used is mainly based on a Luenberger observer. The conversion system is 
divided into two parts, the first mechanical part containing the turbine and the 
gearbox, the second electrical one consisting of a double fed induction 
generator, linked on the stator side directly to the grid, and on the rotor, side 
linked to the grid through two power electronics converters interposed with a 
direct current (DC) link. The machine-side converter is used to control the active 
and reactive powers, and the second on the grid side is used to control the DC 
link voltage. The converters are controlled by the sliding mode strategy, and the 
validity of the methods is checked by simulation using MATLAB/Simulink. 
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NOMENCLATURE 
les Irsc: are respectively the current at the converter grd, prq: _ the rotor d and q axis fluxes. 
outputs on the grid side, rotor side converter 
Vsd, Vsq : the dq axis stator voltages. Rs, Rr: stator and rotor resistances. 
Isd, Isq : the dq axis stator current. Ss, oor : the supply and rotor angular frequency. 
psd, psq : the Stator d and q axis fluxes. E: is the DC voltage that comes from the DC link. 
Vrd, Vrq : the dq axis rotor voltages. Vis2,3: are the three-phase system of the source (the grid). 
Ird, Irq : the dq axis rotor current. Vab,c: are the single voltages from the converter. 


1. INTRODUCTION 

In the last decade the world has experienced a significant development in technology, the latter being 
accompanied by a considerable need for electrical energy. Known as clean energies and given the climate 
changes that the earth has been experiencing since the 19th century, renewable energies are the most used to 
produce electrical energy. This justifies the number of works that have been carried out with the aim of 
improving the quality of the produced energy. 
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are the line currents coming from the source. 
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Among these energy sources, we cite wind energy which has been widely used thanks to the 
advantages it presents. Different structures have been proposed to convert wind energy, there are structures 
based on an asynchronous machine with a permanent magnet, asynchronous machine, and double fed induction 
generator [1]. In the present work, the double-feed induction generator structure is used. This doubly fed 
induction generator (DFIG) generates active power in accordance with the algorithm for maximum power point 
tracking which is widely developed in [2]-[5]. Due to controlling produced power, different control laws have 
been developed in [6]-[15], in our case, we will be interested to the adapted control based on sliding mode thanks 
to the advantages it gives [16]-[20]. 

The major advantage of the wind turbine lies in its resistance to disturbances that can harm its 
operation. we quote [21] who were interested to developed a method to diagnose potential electrical faults 
attacking the DFIG, in [22] the autors focused on the development of a robust algorithm using the sliding mode 
control to avoid the problem of chattering. To reduce this phenomenon the autors of [23] have combined the 
sliding mode and the backstepping controller. For an electrical grid in presence of uncertainty the Hassan et al. [24] 
tried to extract the maximum energy and the control of active and reactive power exchanged. The continuous 
need for energy requires the wind turbine to be put into service for a long time, this causes heating at the level 
of the rotor windings which results in a variation that reaches up to 70% of the nominal rotor resistance, this 
influences produced power which produces less than expected. Therefore, a command adaptation is necessary 
to remedy this problem. 

A set of works have been carried out in the literature to identify the DFIG's parameters, we cite [25] 
who was interested to the rotor and stator inductors identification of DFIG, in Belmokhtar et al. [26] gave all 
the formulas for calculating the stator and rotor resistance. The added value of Boulmane ef al. [27], results in 
the estimation of the stator resistance of an induction machine to reduce the error estimation of its rotational 
speed. Indeed, different techniques have been proposed to estimate an electrical quantity [28], we cite the 
Kalman filters, the Luenberguer observer, and other methods based on the injection of low or high-frequency 
currents (known as the signal processing method). In this work, an adaptation of the control law based on the 
sliding mode technique and using a Luenberger observer is established. To have a readjusted rotor resistance 
estimeted value, a readjustment mechanism using the genetic algorithm is developed [29]. 

This paper is structured as follow. First section is devoted to the studied system modeling with its both 
parts, mechanical and electrical one. In the same section a control law is drawn up based on the sliding mode 
and its robustness is tested. Second section is entirely devoted to the presentation of the DFIG adapted control 
during its parameter’s variation, Luenberger obsever is used to identify this variation. In the first and second 
sections a common part is devoted to the presentation of the procedure used for the optimization of the two 
control approaches using the genetic algorithm. The last section of this paper is devoted to the presentation and 
analysis of the simulation results using the MATLAB/Simulink environment. 


2. RESEARCH METHOD 
2.1. Modelling and control 
2.1.1. General structure of the studied system 

The wind energy conversion systems’ basic structure is given in Figure |. It consists of a DFIG linked 
directly to the grid from its stator when to the rotor, it is through two power converters (C1 and C2) interposed 
with a direct current (DC) bus, and a filter modeling the grid line as seen in Figure 1. The converter switches 
use pulse width modulation of the reference voltages for their control. The reference voltages come from the 
controller blocks, the first block is used to regulate the power produced at its setpoint, and the second controls 
the DC link voltage and imposes a unity power factor. 


Figure 1. Schematic diagram of wind energy conversion system (WECS) based on DFIG 
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2.1.2. Wind turbine model 
The aerodynamic power applied to the turbine blades is given by (1). 


1. 
Poer = 5 Ppsv® 


Where V is the wind speed, p is the air density and S represents the swept area. 
According to betz’s theory the captured power by turbine Pt is given by: 


Pr = Poer Cy(B,A) 
with C, is a specific coefficient to each wind turbine. The turbine torque is defined by: 


P 
Tr = 
Or 


the gearbox is a speed adapter from that of the turbine to that of the generator. Its gain is given by: 


-_ Tr us Omec 


Tmec Or 


Finally applying the fundamental relation of the dynamics, the model is completed by (5). 


dQ).mec 
i 7 ee = Tec — Tem — Tris 


2.1.3. The DFIG model 


() 


(2) 


(3) 


(4) 


(5) 


The double fed induction generator is a very complex non-linear system. To develop a high- 
performance control of this machine [10], [11], for its different operating regimes, precise mathematical 
modeling is required. It can be modeled by its state representation or by equation using the Park transformation. 


The DFIG model in Park's frame is given by: 


Vea = Relea + 24 — W5G5q 
Voq = Relsqg + 22 + ws sa 
Vra = Rylra + 24 — 0, Qrq 
Vig = Rely +22 + 0, Pra 
: = Vealsa + Veglsq 
Qs = Veqlsa — Vsalsq 


2.1.4. The DFIG control 


(6) 


(7) 


To control the DFIG, decoupling is required to facilitate the development of the command, hence the 
need to use the machine stator flux orientation, the latter is oriented along the axis d while cancelling its 
component along the axis q, this is justified by the assumption of the choice of a grad of constant voltage and 
frequency [30], [31]. Taking this hypothesis into account and neglecting the stator resistance (since the DFIG 


is of great power) the equations given previously become: 
Va =0 
Voq = Vy = WsQs 
dor, 
Via = Rrlra + and — Or Prq 
dor 
Veg = Ry lrg ag ae + 0-Pra 
Psa = Ps = Lslsq + MI ra 
0 = Lgleqg + Mg 


Pra = bylra + M154 
Prq = Lrg + MI <q 
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lea = F2 Tl 
ly Ls 
I =] (10) 
sq Te rq 
P; = Veqlsq L Voleg 
(11) 
Vs~s MV 
Q; = Veqlsa = _ = Ty, ira 
M?, dl, M2 
Y= Bane # (eM et gy (br) ig 
2. aly ; (12) 
M dl, M2 MV; 
Vag = Rrleg + lr — To) Ge + 9s (1, - =) Ira + 9 


The control of the DFIG is done indirectly through the converter connected to its rotor (C,). From 
(11) the control of the active and reactive power of the machine is done through the rotor currents by acting on 
the rotor voltages, as illustrated in Figure 2. 


Figure 2. The rotor side converter C1 control 


The basic idea of sliding mode control is to force and pull the dynamics (state) of the system to a 
suitably selected region called the sliding surface, and then to design a control law that will always keep the 
system in this region [22], [32]. From the (11) we deduce as (13). 


Ls Vs (13) 


Knowing that the powers must follow their references given by P™®S and Q’®! we deduce the currents 
reference. 


s(P) = Ie! — Inq 
ref (14) 
s(Q) = La ™ trd 
The control vector is given by (15). 
Veg - V-deq + Veacom (15) 
Vea —_ Vedeq a Vedeom 
To calculate this vector, the derivative of the selected sliding surface is used, which gives as (16). 
: Ls 3 1 MVs 
S(P) = — FE PPh — (Veg — Rrlrg — gWsOLrlra — 9) = —vysgn(S(P)) _ 
1 


S(Q) =- ae + a a i a Rylra + gw,oLylrq) = —v2sgn(S(Q)) 

The choice of the sign function for the sliding surface derivative is justified by the attractiveness 
condition of Lyaponov which requires that the product of the surface and its derivative must be negative or 
zero. Therefore, the constants v,and v2must be positive, as a tool to fix their value, the genetic algorithm was 
used, to have the smallest possible deviation between the power produced and its reference. After developing 
the (16), the control vector is given by (17). 
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Vig = — “ee + Relpg + GWslrOlra + gs + L,ov,sgn(s(P)) 


oe (17) 
V-q = Lo (- ae + ae) + Relrg — gwslyOlyg + L,ov,sgn(s(Q)) 


2.1.5. Parameter’s optimisation using genetic algorithm (GA) for sliding mode technique 

As a method for determining the parameters of the sliding mode controller, the GA was used to have 
a good follow-up of the produced power [33]. In fact, the GA is an optimization and search tool that is based 
on mechanisms that are like the natural selection found in genetics. The GA theory emulates the natural process 
where the most skilful individuals in a population are winners and reproduce themselves while the weakest 
individuals disappear. The method consists of an initial population that is formed by a set of possible solutions 
for the problem, which are coded as chromosomes. These initial chromosomes are submitted to three genetic 
operators: selection, crossover, and mutation. Each chromosome represents a possible solution of the problem, 
and each bit or set of bits, represents a value associated to some variables of the problem. The solutions are 
classified by a fitness function that plays the role of the environment in the natural process [27], [34]. 

As mentioned before, the control system must be able to guarantee the two main objectives: closed 
loop stability of the system and zero regulation of the tracking error. The problem which arises at this level is 
the determination of the sliding mode adjustment coefficients v, and vz. The performance of the control system 
depends essentially on the regulation parameters which are, in our case, the coefficients of the sliding mode 
regulators v,and v2 that we will determine by genetic algorithm. This will make it possible to have a stable 
closed loop system and a good reference's values tracking. The constants v,and v, are fixed using the genetic 
algorithm. They are chosen in such a way as to minimize (12). Indeed, the fitness function (the input argument 
of the main GA program) is defined by the error presented by the sliding surface which must be zero. A GA 
based sliding mode control optimization can be represented as follow Figure 3. The GA properties are given 
in Table 1. The Parameters’ evolution trought the generation is presented in Figure 4. The first simulation 
results are intended to validate the developed command, the simulation is tested for a 1.5MW turbine with a 
variable wind profile. 


Figure 3. GA based sliding mode optimization 


Table 1. Genetic algorithm parameters 


Property Value 
Variables to optimize 2 
Population size 50 
Maximum size of generations 500 
Mutation fraction 0.01 
Crossover fraction 0.08 
Tolerance 5.10? 


Variables’ values 


i i 
oO 10 20 30 40 50 60 7o 80 
Population 


Figure 4. Parameters’ evolution through the generation 


Indonesian J Elec Eng & Comp Sci, Vol. 29, No. 1, January 2023: 120-132 


Indonesian J Elec Eng & Comp Sci ISSN: 2502-4752 m) 125 


Figure 5 illustrates the evolution of the produced active power, which follows the same variations as 
those of the wind profile. It is clearly remarkable that the produced power faithfully follows its reference in 
the case of a DFIG with constant internal parameters. On the other hand in Figure 6, with a rotor resistance 
reaching 50% more than its nominal value, the power produced no longer reaches maximum values [35]. 


108 Active Power 
T 


0 0.5 1 1.5 2 25 3 3.5 4 4.5 5 
Time(s) 


Figure 5. Produced active power with R, = Rn 


x 10° Active Power 


0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 
Time(s) 


Figure 6. Produced active power with R, = 1.5 * Rp» 


3. ADAPTED CONTROL USING THE LUENBERGER OBSERVER 

The full order Luenberger observer is a deterministic state of observation estimating stator currents 
and rotor fluxes. In other words, it is an expression of the model of the system in the state representation which 
makes it possible to reconstruct the evolution of their internal states from measurements of the entry and exit 
of the real system [36], [37]. The basic structure of the control using the Luenberger observer is shown in 
Figure 7. 


Observer Bloc 


Figure 7. The RSC Adapted control diagram using a Luenberger observer 
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3.1. Dynamic model of the DFIG in the coordinate (a,B) 

To estimate the value of the rotor resistance which changes during the DFIG operation, a state 
representation of the latter must be made. It can be described by the following state equations in the stationary 
reference (a, B): 


= [A][X] + [B][U] 
[ 


18 
Y] = [C][X] (18) 
Ise 
I 
1 =| 9 (19) 
drp 
Vsa 
ei | (20) 
= 0 6 0 
[B] = es & (21) 
ote B 
Isc 
IY] = is (22) 
—A Ws £ Bo 
_|-e@s -A Bo f 
lal=| : = 
as 0 =f -w 
TR TR 
& 2 4 & 
TR TR 
1000 
=|, 4.6 (24) 
M L M2 
Wit! a oO Taha 


According to Kalman a system is observable if the rank of the matrix y is equal to n with n being the 
dimension of the state vector of the system, in this case n = 4. The matrix y is defined by (25). 


C 


_1CA 


CA? 
After calculation we find: rang[w] = 4. This result confirms the observability of the DFIG. 


3.2. Luenberger observer structur 
The Luenberger observer structure is illustrated in Figure 8. The Luenberger observer equation is 


given by (26). 
[| = (4][X] + (Bl(u] + (LILY - 7] 


(26) 
[7] = [cllx 


The estimated quantities from this observer will be used to reconstruct the value of the rotor resistance, but 
long before it is necessary to determine the elements of the matrix [L] defined by (27). 


Ly —L 
_|l2 la 

=|? 7, (27) 
ly Ls 
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Figure 8. Luenberger observer structur 


To do this we define the following matrix: 
[4] = [4] - [z][c] (28) 


the technique used is that of pole placement, we start by solving the characteristic equations of the full-order 
observer and the system’s one expressed by: 


det(SI — A) =0 (29) 


with J identity matrix. After solving the two equations the gain matrix L is calculated so that the observer poles 
are proportional to those of the DFIG. After calculation and identification, we find (30). 


L,=-(k-DA+t 


Ly =k-1 
= Rt(, _ MB) k-tey y 4 (30) 
Ls =" ( ve) + Sats) 

ee 

Ly =~ 


3.3. The adaptation mechanism of Rp calculation 
Initially, the quantity to be measured must appear in the state matrix of the system and the observer, 


so the reformulated state matrixes are given by (31), (32). 


-—1 a B ~ 
“ig Ws ig RR Bo 
ote —-Ws a Re —Bw FR 
[A](Re) =| wo ais G1) 
Tg RR 0 ink —-W 
M x -1¢4 
0 in RR @ rs 
-1 B 
Sig ARR 0 1g ARR 0 
0 AR, 0 AARg 
AA=| y ae (32) 
—ARp 0 —ARp 0 
LR LR 
M -1 
0 7, ORR 0 7, ARR 


The adaptation mechanism will be deduced by Lyapunov's theory. The stator current estimation error 
and rotor flux, which is the difference between the observer and the motor model, is given by: 
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we consider the following Lyapunov candidate function: 


v = ete + Rk (35) 
RR 


its derivative is given by: 


dv _ det tde, 2 dAR; 


mr dee! ae ie HE (36) 
@ = (de + AAR) ¢e + et(Ae + AAR) — 5a, AR ae (37) 
— = [e'(A‘ + A)e] + [(AAX)‘e + e'AAX] - [ake tie (38) 
the development of the second term is given by: 
(AAR)‘e + e'AAR = 2e'AAK = 2 « |(—ARgheg + -ARe Pra) Cisa + (—ARelsp + 
Zig hse © 7, RP Re) Pigg Lo nse 
B z M A 1 z M A 1 z 
F ARePep) Cig + ( ARalsa = 7 ARR Pra) ebea + ie ARplsg — AR Pep epep| (39) 


The derivative of the Lyapunov candidate function must be negative to check the stability of the 
system, the first and third terms being negative, it remains to verify the second or else impose it zero, neglecting 
the error of the rotor flux we find the following (40) and (41): 


2+ [(St tae + E bac) tg + (ZElbn + Ea) eg] BRe ~72-aRy ME = 0 a) 


oLR 


x = % x a a 
Rr = Opp [FA Sdsa Cisq + Ispeig, dt + £ bra Clgg + Prpeig, At] (41) 


We pose: Cy = Isqei,, + IsBeisp and Cy = Pralicsg + Prpeisg 

To have a stable system instead of estimating the value of the resistance using an integrator, we opt 
for the solution of a proportional integrator. Finally, the rotor resistance adaptation mechanism is illustrated in 
Figure 9. 


Figure 9. Adaptation mechanism structur 


The integrator and the proportional gain of the adaptation mechanism are calculated using the genetic 
algorithm cited before, to have the value of the rotor resistance estimated as close as possible to the real value. 
These parameters are adjusted to obtain the optimum performance in terms of dynamics, robustness, and 
disturbance’s rejection. The objective function is therefore chosen to minimize the difference between the rotor 
resistance mestimated and the reference one. A GA based Adaptation mechanism optimization can be 
represented as follow Figure 10. 
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Figure 10. GA-based Luenberger observer optimization 


4. SIMULATION RESULTS 


Simulations were carried out using MATLAB/Simulink. 1.5 MW wind turbine has been used, its 
parameters are given in the appendix Table 2. During the simulation, rotor resistance was increased by 50% of 


its nominal value. Simulation results for 690 V of stator voltage and DC link reference equal to 1200 V are 
given by these Figures 11-13. 


Table 2. Doubly fed induction generator parameters 


Parameter Symbol Value 
Rated power P, 300 Kw 
Frequency f 50 Hz 
Stator resistor Rs 8.9 mQ 
Rotor resistor Rrn 9,13 mQ 
Stator inductor L, 12.9 mH 
Stator inductor L, 12.7 mH 
Mutual inductor M 12.672 mH 
Dispersion coefficient o 0.0198 
Number of pole pairs P 2 


Simulation results compared to the first simulations given in the modelling part, show the monitoring 
of the reference power even in the presence of the variation of the rotor resistance as shown in Figure 11. In 
fact, the obsever identify the rotor resistance value see Figure 12. In Figure 11 that the DFIG produced power 
follows its reference from 300 ms using the control adapted by sliding mode. This delay is because of the 
estimated rotor resistance becomes equal to its reference from the same instant as illustrated in Figure 12. It 
should be noted that the reactive power also see Figure 13 remains zero equal to its reference from 300 ms. 
Before this moment active power, reactive power, present a peak for a short time. 


ox 0° Active Power 
if T T 


0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 
Time(s) 


Figure 11. Active power with rotor resistance variation and estimation 
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Figure 12. Reactive power with rotor resistance variation and estimation 
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Figure 13. Reactive power with rotor resistance variation and estimation 


5. CONCLUSION 

In this work, the study, modelling, and control of a wind energy conversion system based on a DFIG 
were investigated, the evaluation of the dynamic performance of this system controlled by the sliding mode 
strategy and its sensitivity to the variation of the internal parameters of the DFIG. The control approaches used 
to control the power exchanged with the grid consist of adapting the DFIG rotation speed to the wind speed by 
controlling the rotor currents and the electromagnetic torque. We have therefore developed a robust and 
adaptive control strategy based on the sliding mode and the Luenberger observer with optimization by genetic 
algorithm. The latter allows to have optimal values of the controller parameters to guarantee the best 
performances. For normal operation, the internal parameters of the DFIG are assumed to be fixed. From the 
simulation results we note that the sliding mode control technique, optimised by GA, provides good 
performance, including fast response and no overshoot. In practice, the internal parameters are always subject 
to variation due to several physical phenomena. We took as an example the variation of the rotor resistance 
and found, from the simulation results, that the power output no longer follows its reference for the 
conventional method. However, the use of the adapted control, proposed in this paper, based on the Luenberger 
observer presented more satisfactory results. In terms of perspectives, the developed control strategy can be 
further complemented by the implementation of an observation block taking into consideration all internal 
parameters including both stator and rotor resistances as well as inductances. 
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